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The origin of [Fell] emission in NGC 4151 
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ABSTRACT 

The centre of NGC 4151 has been observed in the J-band with the SMIRFS integral 
field unit (IFU) on the UK Infrared Telescope. A map of [Fe II] emission is derived, and 
compared with the distributions of the optical narrow line region and radio jet. We 
conclude that, because the [Fell] emission is associated more closely with the visible 
narrow-line region than with the radio jet, it arises mainly through photoionization 
of gas by collimated X-rays from the Seyfert nucleus. The velocity field and strength 
with respect to Pa/3 are found to be consistent with this argument. The performance 
of the IFU is considered briefly, and techniques for observation and data analysis are 
discussed. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: jets - galaxies: nuclei - 
instrumentation: spectrographs - techniques: spectroscopic 



1 INTRODUCTION 

Seyfert galaxies exhibit strong emission lines from forbid- 
den electronic transitions — notably thos e of oxygen, n itro- 
gen and sulphur at visibl e wavelengths (Seyfert 1943) and 



iron in the near-infrared (Rieke & Lebofsky 1981; Simpson 
|et al. 1996 ). These transitions can occur where ions are ex- 



cited by collisions with electrons, yet can decay r adiatively 
before subsequent colli sions lead to de-excitation (Peterson 



low-der 



1997; QsLurbroek 1989). Hence the lines are associated wlLh 



fclly it 



regions ufiuiiUud hydiugeii, wheie lieu electiuiis 

are abr Jndant but the timescale for collisions is suffici e ntly 
long. 

Some species, including neutral oxygen and Fe + , have 
similar ionization potentials to hydrogen, so do not survive 
where it is fully ionized. The most favourable conditions for 
[OI] and [Fell] emission arise in partially ionized zones, the 
extent of which influences the observed intensity. Since H II 



regions around stars normally give way sharply to neutral 
gas beyond the Stromgren radius, their [O I] and [Fe II] lines 
are relatively weak. The central spectra of Seyferts, in con- 
trast, indicate the presence of extensive partially-ionized hy- 
drogen clouds as well as very highly-ionized species. These 
differences are explained somewhat naturally if the emis- 
sion line clouds in Seyferts are photoionized by the 'power- 
law' continuum of the central active galactic nucl eus (AGN), 
whic h is spectrally flat ter than a black body ( Osterbrock 
198S i_ Mouri et al. 1990|). Numer ical simulations ( Ferland & 



Netzer 1983; Mouri et al. 2000) lend credence to this see- 



The second mechanism which can generate enhanced 
zones of [O I] and [Fell] emission is shock excitation of gas 
( |Osterbrock 1989 ). This appears to be the dominant cause in 
starburst galaxies, where shocks due to supernova remnants 
are important and sup ernovae have been associated directly 
with enhanced [Fell] (Lester et al. 1 99Cj; Greenhouse et al 



Forbes et al. 1995 



Greenhouse et al. 1997). Compar- 
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1991 

ison of observed spectral line ratios with recent numerical 
models helps to distinguish between shocking and photoion- 
ization (Mouri et al. 2000), and indeed proves consistent 
with the former dominating in starbursts and the latter in 
Seyfert nuclei. Doppler broadening can also be used to check 
for kinematic disturbances. 

Long-slit spectroscopy of Seyfert galaxies has revealed a 
spatially extended narrow line region (ENLR) in the visible 
(at radius greater than 4 arcsec in the case of NGC 4151) 
beyond the central narrow-line region (NLR). Unger et al. 
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(1987) traced the bright [OIII] lines at 4960A and 5007A, 
generated in fully-ionized regions, along with H/3 (486lA) in 
a sample of galaxies. The features were shown to have mini- 
mal lin e-of-sight velocity structure consistent with simple ro- 



tation, ind line ratios characteristic of high excitation levels, 



poi nting once again to photoionization. Subs equent stud- 



ies (Penston et al. 1990; Robinson et al. 1994) have tended 



to support this conclusion. The elongated, sometimes visi- 
bly conical morphology of ENLR (and NLR resolved with 
the Hubble Space Telescope: HST) poi nts to a link with 
collimated radiation from the nucleus ( Unger et al. 1987 



Mulchaey et al. 1996; Evans et al. 199; ) 



In the near-infrared, the origin of strong, extended 
[Fell] emission in Seyfert galaxies has been somewhat con- 
troversial, with evidence for both shock excit ation and pho- 
toionization occurring. It has been argued (Greenhouse et 
al. 199lj) that shocks must be important because iron in our 
interstellar medium is mainly co ndensed onto grains whic h 
shocks can destroy via sputtering (Savage & Sembach 1996) 



However, recent studies (Simpson et al. 1996 



Mouri et al 



200C|)jJidicate that the ratio [FeIIJ/[OI] in both Seyferts 
and starbursts is consistent with the usual depletion levels 
of iron and oxygen from the gaseous phase, and that certain 
other signs of grain destruction are absent. Even if metal 
enhancement is not the dominant factor, it is likely that 
shocks due to outflows along the radio jets of AGN do make 
s ome contribution to emission from p artially-ionized zones 
( Simpson et al. 1996[ [Morse et al. 1996| ). In particular, Knop 
et al. 1996 find that whilst [Fell] in NGC4151 seems kine- 
matically compatible with the ENLR, the 1.257 (jm feature 
is broader than the nearby Pa/3 line 

NGC 4151 is both a natural candidate in which to study 
[Fe II] excitation and an important test case for understand- 
ing the structure of active galaxies. It is one of the brightest 
nearby Seyfert galaxies (V ~ 11.5, cz — 995 kms" 1 ), yet 
eludes firm classification in the framework of unified AGN 



models (Antonucci 199E ) . The highly extended narrow line 



emission and radio jet ( Perez et al. 1989 ; Pedlar et al. 199; ) 



are characteristic of a Seyfert 2 nucleus, whose collimation 
axis is almost perpendicular to the line of sight. However, 
the object also exhibits the strong broad lines and variability 



of a Seyfert 1 (Maoz et al. 1991), suggesting quite a differ- 



ent orientation. NGC 4151 has therefore been the subject of 
many past projects involving slit spectroscopy and narrow- 
band imaging. 

An unusual feature which facilitates the present inves- 



at position angle ~50° ( 


Perez et al. 1989; Evans et al. 1993 


Kaiser ct al. 200C), and radio iet at 77° ( 


Pedlar et al. 1993 



Mundell et al. 1995 ). Whatever its origin, this provides an 
opportunity to associate [Fell] directly with the optical 
ENLR and/or radio jet, by tracing the emission along both 
axes. If a clear separation can be made, it will be possible to 
determine the relative influence of ionization by collimated 
X-rays and shocking by outflowing radio plasma. Such a 
task is inherently suited to integral field spectroscopy, which 
can form a complete, homogeneous picture without prior as- 
sumptions, slit alignment problems or the waveband limita- 
tions of a filter. 

The intrumentation is discussed in the next section, fol- 
lowed, in Section 3, by details of the observations. In Sec- 
tion 4, we present the method of data reduction and analysis 



which presents particular challenges since the technological 
limitations required the fibre outputs to overlap at the slit. 
However this does not degrade the data provided that the 



field is critically sampled at the IFU input ( Allington-Smitb 



& Content 1995 ). Finally, the results are discussed in Section 



5 and compared with other work. 



2 INSTRUMENTATION FOR INTEGRAL 
FIELD SPECTROSCOPY 

The S MIRFS-IFU is an ex p erimental 72-fibre integral field 
unit ( |Haynes et al. 1999|; |Lee 1998|) feeding the CGS4 
spectrograph Wright et al. 1993] ) at the UK Infrared Tele- 
scope (UKIRT). It provides integral field spectroscopy in 
the J and H bands via lensed fused-silica fibres. It was 
designed and built by the Astronomical Instrumentation 
Group at Durham University as a prototype for much larger 
integr al field units: TEIF U on the William Herschel Tele- 
s cope(|Murray et al. 200o|) a nd the GMOS-IFU on Gemini 
( Auington-Smith et al. 2000). Futher details may be found 
in Haynes et al. (1999). 

SMIRFS is designed to intercept the telescope beam 
without modification to existing instrumentation, but con- 
sequently requires removal of the CGS4 calibration unit. In 
place of this module, the SMIRFS slit projection unit feeds 
light from the fibre bundle into the long-slit, circumventing 
the usual beam path; the IFU input is connected to another 
port of the UKIRT instrument support structure. As a re- 
sult of this arrangement, any reference observations must be 
taken using the main telescope aperture. 

At the image plane, a hexagonal microlens array pro- 
vides optimal coupling to the slow telescope beam and avoids 
dead space between fibres (see Haynes et al. for more infor- 
mation). The field of view is 6 x 4 arcsec, with 0.62 arcsec 
spacing between adjacent fibre centres. A second set of mi- 
crolenses at the output restore the correct focal ratio for 
CGS4. When using the long camera, each fibre projects a 
spot whose FWHM is ~1.8 pixels at the detector, giving a 
comparable spectral resolution to the standard 2 pixel wide 
slit. These spots are spaced two pixels apart, placing the in- 
strument on the borderline between the 'maximally packed' 
and 'resolved peaks' categories defined by Turner (2001), al- 
though throughput variations obscure the regular pattern 
in practice. Since there is only one field, sky subtraction is 
performed by nodding the telescope off source. 

The capability of the SMIRFS-IFU was previously 
demonstrated in a June 1997 commissioning run. Its 
throughput relative to CGS4 alone was measured to be 
~50%; since then, minor adjustments have been made to 
improve performance. 



3 OBSERVATIONS 
3.1 Science programme 

NGC 4151 was observed with the SMIRFS-IFU on 15- 
17 February 1998, as part of a sample of Seyfert galax- 
ies. Details of the measurement parameters are presented 
The selected region of the J-band includes 



in table 3.1 



both [Fell] and Pa/3, at rest wavelengths of 1.2567^im 
and 1.2818/xm respectively. Previous studies have sometimes 



© RAS, MNRAS 000, ??-?? 



inNGC4151 3 



Table 1. Details of NGC4151 observations 



Instrument 


CGS4 + SMIRFS-IFU 


Camera 


300mm 


Grating 


150 lines/mm 


Diffraction order 


3 


Filter 


Bl 


Spectral range 


1.25-1.31/Ltm 


Spectral FWHM 


~ 4.2A (lOOkms -1 ) 


Spectral sampling 


2.lA 


IFU field 


6 x 4" 


Spatial FWHM 


~ 1.2 x 1.1" 


Spatial sampling 


0.62 x 0.54" 


Detector 


256 x 256 InSb 


Gain 


6e — /data number 


Read noise 


23 e~ 


Dark current 


~ 1 e"/s 


Single exposure 


180s X 2 detector positions 


Central exp. 


48 minutes on source 



used [Fell] at 1.6435/im and Br7 at 2.1655/xm; such mea- 
surements are equivalent to ours, because the same upper 
energy level of Fe + is involved. Since, however, the longer- 
wavelength lines are much further apart, they cannot be 
observed together at high dispersion and their ratio is more 
sensitive to differential extinction along the line of sight. 
The J-band [Fell] feature is also brighter, and the uncooled 
SMIRFS better suited to working in that regime. 

The spatial resolution quoted in table 3.1 is for the fi- 



nal, combined data set; it was estimated by producing im- 
ages of the unresolved [S IX] line and Pa/3 broad line. These 
FWHM values correspond to the Nyquist scale, suggesting 
that the seeing disc was undersampled — the implications are 
discussed briefly in 

The long axis of the IFU was orientated at 50° east of 
north, along the optical ENLR of the galaxy. The field of 
view, h owev er, only covers the inner NLR at < 4" (Wingc 
et al. j999|), which is orientated less favourably at ~60° 



(Kaiser et al. 2000). With arcsecond resolution, this gives 
~1 FWHM separation from the radio axis. Hence a mosaic 
was made of several telescope pointings around the nucleus. 

Target centering was performed with the aid of a dedi- 
cated program, which derives a rough map of the illuminated 
input elements from the data. For the peripheral mosaic po- 
sitions, it was necessary to rely on specified telescope offsets 
from the centre. Each position was observed in turn, fol- 
lowing the repeated sequence object-sky-object; this maxi- 
mizes the time on source without giving unequal object and 
sky exposures at individual pointings. The 180s integration 
time was just long enough for the background to dominate 
the read noise and dark current. The whole mosaic was re- 
peated until the full exposure was reached. 

In the absence of the calibration unit, sky emission lines 
were used as a wavelength reference. Some arc lamp obser- 
vations were also taken through the telescope. The IFU was 
flat-fielded using the illuminated dome and the twilight sky, 
whilst a detector flat was made with the calibration unit, 
prior to the installation of SMIRFS. The bright F and A 
stars BS4572 and HD105601 were observed in order to cal- 
ibrate the spectral response of the instrument and measure 
sky absorption features. 




Figure 1. Sketch of the SMIRFS-IFU input with a mask in place 
for determining the output magnification. The line with an arrow 
indicates the fibre ordering at the slit. The dead fibre is shaded. 



3.2 IFU set-up and characteristics 

The IFU was aligned at the focal plane by projecting a laser 
beam back through the fibres onto the UKIRT secondary 
mirror. The CGS4 slit was set parallel to the row of fibre 
outputs, before adjusting the optics of the slit projection 
unit so that fibres were imaged two pixels apart in the spec- 
trograph. The two pixel magnification was verified by illu- 
minating the IFU with a mask in place, generating an image 
of the fibre slit at the detector with a gap for every input 
row. Since the mapping of the input to the slit reverses every 
row, to ensure that elements are adjacent both at the input 
and output (Fig. |l|), the distance between illuminated band 
centres at the detector alternates about the row separation. 
By fitting a straight line to the centroids, the magnification 
was measured as 2.002 ± 0.005 pixels/fibre. 

The J-band throughput of the IFU was measured by 
comparing observations of the standard stars HD 84800 and 
BS 4069, taken with the CGS4 slit and the SMIRFS-IFU re- 
spectively. Examination of the point-spread function (PSF) 
indicates that truncation of the seeing disc by the 4" slit 
did not cause significant light loss. The difference in airmass 
was small, giving an expected error due to atmospheric ex- 
tinction of < 1%. It was found that the throughput is 50%, 
consiste nt with th e figure of 49% from the 1997 commission- 
ing run ( Lee 1998 ). 

Variations in the fibre transmissions were quantifie d by 
extracting a flat-field spectrum for each element (see §Kk2) 
and integrating in wavelength (Fig. ^). The values there- 
fore include crosstalk from neighbouring fibres. For Gaussian 
profiles with FWHM 1.8 pixels, the contribution from either 
neighbour is ~10% of the total. Ignoring further smoothing 
due to interpolation and any correlation between through- 
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Element number 
Figure 2. Throughput variation along the slit. 



puts, the measured and true RMS values are thus related as 
follows: 



""data 



a/(0.8cti FU ) 2 + 2 x (O.Icttfu) 2 = 0.8cri FU . 



The measured RMS, <7data, indicates the signal-to-noise fluc- 
tuation; it is 8% overall or 6% excluding the dead fibre and 
some vignetted elements. The corresponding values of crnro, 
which characterizes the fibre bundle, are ~10% and ~7% 
respectively. It is concluded that whilst careful flat-fielding 
is important, inherent differences in the signal-to-noise ratio 
are reasonably modest. 

The point spread function does not var y substantially 



bet ween fibres. Whilst a previous estimate (Haynesetal 
199s()~i}idicat es that some profiles are up to 20% broader 
than others, the PSF is only known to an accuracy of ~1Q%. 
Any variations at this level will contribute slightly to the 
Nyquist-scale noise in line width maps. The effect on line 
flux noise is negligible, since the broadening is a fraction of 
a pixel. Any inhomogeneities introduced by the fibres can 
be smoothed out to some extent by combining observations 
with slight spatial offsets. 



4 DATA REDUCTION AND ANALYSIS 
4.1 Overview 

Data were reduced in IRAF, using several purpose-written 
scripts and programs alongside standard tasks. As discussed 
in Turner (2001), the best method of dealing with fibre 
spectra depends largely on their separation and FWHM at 
the detector. For the SMIRFS-IFU, spectra are much more 
densely packed than in a multi-object design, precluding 
the use of existing extraction software. New tasks were also 
needed for reconstructing and working with the observed 
x-y-X volumes in the form of datacubes, and for bad pixel 
correction. Because the slit maps to the field of view in a 
complex way, flexure in the spectrograph cannot be compen- 
sated by re-centring the target (this would be true anyway 



for off-centre mosaic positions). Observations must there- 
fore be co-added after conversion to datacubes, rather than 
in their raw format. This requires effective bad pixel correc- 
tion for integration pairs. 

Turner (2001) presents an alternative reduction method 
for fibre-IFU data, inspired by the very dense pac king of 
spectra, one per row, in TEIFU (Murray et al. 2000). With 



the long CGS4 camera, the raw format of SMIRFS-IFU data 
is intermediate between that of TEIFU and multi-object de- 
signs. Since fibres are imaged with a FWHM of 1.8 pixels, the 
output pattern is almost critically sampled and is acceptable 
to interpolate in pixel coordinates. Moreover, the two pixel 
fibre separation is large enough that moderate smoothing by 
the interpolant will be small in relation to the image reso- 
lution. At the same time, flexure-related errors and under- 
sampling at the IFU input rule out much gain in precision 
through TEIFU-style reduction. Hence, as for most fibre- 
based instruments, the reduction is based on extracting ID 
spectra. 

Finding absolute fibre positions on the detector after 
flexure is complicated by vignetting at the slit ends, which 
prevents direct location of the first and last elements. With 
an already limited field of view, it was decided to rely on 
the single broken fibre for measuring offsets between frames, 
rather than masking off row ends at the input to create 
regular gaps. This is usually adequate, if somewhat awkward 
and instrument-specific. 

With 0.62" elements at the IFU input, undersampling 
is potentially an important issue. Inspection of various ob- 
servations (with and without SMIRFS) indicates that the 
seeing during the February 1998 run was ~0.9", including 
integration over the input microlenses; this figure is rel- 
atively favourable for such coarse sampling. A very basic 
simulation of spatial frequency propagation through the in- 
strument shows that, depending on the alignment of fibres, 
pixels and image features, aliased frequenci es in the dat a 
can have amplitudes up to 10% of the mean ( Turner 2001 ). 
However, the Fourier spectrum is only affected significantly 
above ~80% of the Nyquist frequency, so moderate smooth- 
ing during extraction and reconstruction can easily eliminate 
much of the contamination. In principle, offset observations 
can be combined w ith phase shifts such that aliasing is re- 
moved (Lauer 1999), but in this case the main effect would 
be to amplify noise because the offsets are not known or con- 
trolled precisely (and the seeing varies). Nevertheless, simply 
co-adding datacubes after alignment on a common grid does 
suppress residual alias components, which add out of phase 
like noise. At the end of the reduction process, aliasing ar- 
tifacts are not a major source of error. 

4.2 Calibration and extraction 

Each observation consists of two integrations, offset by 
one detector pixel in the spectral direction to account for 
bad pixels. These raw frames are largely obscured by non- 
uniform dark current. The first step was to subtract a sky 
pair from an object pair, removing dark current in the pro- 
cess. However, this double use means that sky frames cannot 
be smoothed spatially to reduce their noise level before sub- 
traction. The result was then divided by a detector flat-field 
image, obtained with the calibration unit before the IFU 
was installed. 
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A less accurate dark frame, taken earlier in the night, 
was also subtracted from the object separately. This reveals 
the fibre throughput pattern so that the dead fibre can be 
located. Ideally, this would be done using a sky minus dark 
image, but flexure of a fraction of a pixel can occur even 
between consecutive pointings. This is important for image 
reconstruction, but less so for sky subtraction — assuming 
Gaussian fibre profiles and an offset of 0.3 pixels, the sky 
residual due to throughput variations is likely to be around 
3%, falling below the noise level. 

A program was written to combine the pairs of offset, 
sky subtracted frames with bad pixel removal. This makes 
a crude noise estimate based on the RMS residuals between 
the two images; any input pixel which differs from the me- 
dian of its neighbours by more than a specified number of 
standard deviations is excluded. The procedure can some- 
times truncate narrow arc lamp emission lines, whose gra- 
dients can be much greater than the noise level, but is safe 
for the galaxy spectra (which have had sky lines subtracted 
already). Any obvious remaining blemishes were removed 
manually. 

Spectra were straightened with respect to detector rows 
by applying a small pre-determined shift to each column. A 
curvature model was constructed by measuring the centroid 
of a suitable image at each wavelength and fitting a low- 
order polynomial to the values. Interpolation was performed 
using a cubic spline, which is good for dealing with slightly 
undersampled data. Inspection of corrected frames verifies 
that the curvature is practically constant along the slit and 
hence independent of flexure. The maximum difference be- 
tween columns is 0.7 rows. 

Absolute fibre positions were found by summing the 
straightened object minus dark image in wavelength and fit- 
ting the dead fibre position using standard IRAF routines. 
The corresponding sky subtracted image was resampled so 
that fibres fell at the mid-points between odd pairs of rows. 
These pairs were combined with the blkavg task, producing 
a stack of extracted fibre spectra (boundary rows were dis- 
carded). Sometimes spatial gradients make the dead fibre 
difficult to identify, but the procedure can be repeated if a 
reconstruction appears distorted. 

The stack of extracted spectra was divided by a fibre 
flat, created in the same way but averaged over its wave- 
length range. Sky emission lines in the object minus dark 
image were used for wavelength calibration. These allow 
correction of grating offsets between observations and also 
provide an absolute calibration. The spectra from every ob- 
servation were interpolated onto a suitable standard grid, 
truncating the end few pixels which fell outside the com- 
mon range. The error in determining the relative offsets is 
around 0.2 A. The error in the linear wavelength fit is 0.7A, 
which is more than adequate for the present work. Separate 
solutions were not determined for the spectra within each 
observation; variation along the slit (±0.3A) is smaller than 
the final random velocity errors in the mosaic and much 
smaller than the true gradient in the galaxy. 

Finally, the stack was divided by a continuum- 
normalized stellar spectrum, to remove telluric absorption 
features. A Pa/3 line in one of the stars was interpolated 
over. Unfortunately, the mosaicing process makes it diffi- 
cult to re-observe the standard stars frequently, so there 
was sometimes a non-negligible difference in airmass from 



the galaxy. Although the main absorption band at 1.269/im 
has been removed effectively in the final data, there may be 
some absorption residuals at around 2% of the continuum 
level. 



4.3 Construction of datacubes and mosaicing 

Sets of extracted spectra were interpolated spatially at each 
wavelength step, onto a square grid of 0.15" pixels.^] The 
constituent datacubes of each mosaic cycle were combined 
with appropriate offsets using the program mosaic, which 
has been inclu ded as a prototype task in the imspec package 
(Turner 2001). This calculates the smallest common grid 
from a list of pointings and resamples the images at each 
wavelength using a tapered sine function. The precise inter- 
polation method is not important, since the datacube grids 
are much finer than the resolution of the data. A mask im- 
age is also produced, recording the number of input pixels 
contributing at each point. 

Unfortunately, telescope pointing errors of a fraction of 
a fibre are evident. This is partly because a dichroic mirror 
is used to guide on the target at visible wavelengths whilst 
observing in the infrared; the two images drift apart as atmo- 
spheric dispersion changes with the airmass along the line of 
sight. Modelling this effect may systematically reduce errors, 
but was outside the scope of the current project. There may 
consequently be some image distortion at the points of over- 
lap; mosaicing nevertheless reveals any wider-scale extension 
of spectral features within the noise limit of the data. For 
the NGC 4151 observations, emission turns out to be con- 
fined mainly to the central field; three mosaic positions with 
appreciable signal-to-noise, offset diagonally east and west, 
were finally used for mapping. 

After constructing each separate mosaic, the resulting 
datacubes were combined using comcubes, also included in 
imspec. This is very similar to mosaic, but calculates the 
offsets between datacubes by locating the centroids of the 
intensity peaks after integrating in wavelength. The individ- 
ual overlap masks are also combined to show the total num- 
ber of data points contributing at each output pixel. The 
final datacube is divided by the mask, to scale the pixels 
correctly. 

The spectral response of the instrument was deter- 
mined by fitting piecewise polynomials to the continua of 
the standard star observations, dividing the results by black- 
body curves appropriate for their spectral types. This was 
achieved using the standard, sensfunc and calibrate tasks in 
the IRAF onedspec package. The response varies at the level 
of 5%, and with atmospheric absorption features already re- 
moved it is not strongly dependent on airmass. The accuracy 
to which the curve can be determined is limited mainly by 
the wings of absorption lines in the stars and atmosphere, 
which are difficult to distinguish from modulation in the 
continuum. The final datacube was divided by the response 
function using calibrate directly and the values normalized 
to a unit mean. 



We used the method of Renka and Cline(Rcnka 1983) for 



convenience. 
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Figure 3. (a) 1.25— 1.3/xm image of the nucleus of NGC4151. (b) Paschen fS broad line image. This indicates the spatial resolution. 

4.4 Mapping the flux and velocity field 

A velocity map was obtained from the combined mosaic. 
Spectra with mean fluxes below 2.5% of the central value 
(about 3<t) were not included, since line emission is barely 
detected at the corresponding distance. The NOAO task 
rvidlines was used to measure a combined [Fell] and nar- 
row Pa/3 redshift for each spectrum. This is more reliable for 
low signal-to-noise peripheral points than using [Fell] alone. 
The results show that differences in the velocity fields of the 
two lines are comparable with the measurement errors and 
small compared to the line widths. Hence any deviations in- 
troduced by Pa/3 will not affect line strength measurements 
significantly. The plain text results from rvreidhnes were 
redirected into a file and converted to a 2D image using 
listoim in imspec. 

Line fluxes were measured using the specmap task in 
imspec. This reads the datacube, velocity map and a list 
of wavebands for line and continuum estimates, producing 
a flux map as output. The measurement bands and linear 
continuum estimates can be overplotted on the spectra for 
verification. Over the wavelength range observed, there is lit- 
tle pure continuum to provide a suitable reference for [Fe II] , 
but experimenting with different bands shows that the de- 
rived spatial profiles are robust. 

An excitation map, showing the strength of [Fell] rela- 
tive to Pa/3, was also produced. Since this involves dividing 
one image by the other, the result is sensitive to noise where 
emission is weak. The map was therefore created using line 
images from a spatially smoothed version of the datacube. 
Rough line width estimates were made by examining indi- 
vidual spectra with implot and splot. 
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Figure 4. Central spectrum of NGC4151. 



5 RESULTS AND COMPARISON WITH 
OTHER DATA 

A wavelength-integrated image of NGC 4151 is shown in 
Fig. [](a). Whilst some horizontal broadening is caused by 
the IFU, because spectra overlap, there is also a real NW-SE 
extension of the continuum, following the stellar distribution 
seen in optical images. Fig. ^|(b) is an image of the Pa/3 broad 
line (including the narrow line at ~6% of the total), which 
gives an indication of the spatial resolution. 

A central spectrum is presented in Fig. ^; this includes 
the broad and narrow Pa/3 components, [Fell] and [SIX]. 
There may also be contributions from a weaker 1.2703/im 
[Fell] feature, superimposed on the blue wing of the Pa/3 
broad line, and He I at 1.2791pm, to the left of the Pa/3 
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Figure 5. Combined [Fell] and Pa/3 velocity map. 



narrow line. The small bump to the right of [S IX] is probably 
due to imperfect removal of atmospheric absorption features. 
Th e combined [Fell] and Pa/3 velocity map (Fig. ||) 



traces t jhe rotation of gas about an axis projected at a posi- 
tion angle of ~110°. It is broadly consistent bo th with pre- 
vious slit measurements of the same features ( Knop et al. 
199f ) a id with o ptical spectroscopy of the NLR ( [Vlediavilla 
& Arri bas 1995] ; [Winge et al. 1999| ). Estimates of the line 



widths are complicated by the Pa/3 broad line shape, lack 
of clear continuum and modest signal-to-noise ratio off cen- 
tre. The line profiles have not been modelled with multiple 
components, but their overall deconvolved widths are close 
to the values found by Knop et al, who note that the width 
of [Fell] relative to Pa/3 suggests a contribution from shock 
excitation. At the centre, the FWHM values are approxi- 
mately 400 kms- 1 for [Fell] and 250kms~ 1 for Pa/3. Both 
lines are broader in the brightest parts of the surrounding 
extended emission (see below), reaching ~600kms _1 and 
^400 km s -1 respectively. Additional broadening in the di- 
rection of the radio jet is not evident. 

As an alternative to Knop et al.'s suggestion that the 
greater velocity width of the [Fell] emission is due to a 
contribution from shock excitation, we note that our result 
i s consistent with t hat for a sample of Seyfert 2 galaxies 
( Veilleux et al. 1997). Veilleux et al. cite this as evidence that 



the two lines are emitted in different volumes of gas. We also 
note that Kaiser et al. (2000) found different components 
in the velocity field of NGC4151 revealed by [OIII] emis- 
sion. They suggested that the low velocity dispersion com- 
ponent arises from smooth flow in the gravitational field of 
the galaxy while the high dispersion component results from 
an outflowing wind. Our result could arise if the Pa/3/ [Fell] 



ratio was higher in the low dispersion component than in 
the high dispersion component. 

The spatial distribution of [Fell] is mapped in Fig. 
alongside examples of individual spectra. The image pro- 
file is extended along both axes of the IFU with an in- 
trinsic RMS of 0.3-0.4" horizontally and 0.9-1.0" vertically 
(FWHM 0.8" and 2.2" respectively assuming a gaussian dis- 
tribution). The narrowest point is at the intensity peak, as 
might be expected for a biconical shape. Although these 
scales are comparable to the image resolution, it is clear that 
the [Fell] is aligned closer to the optical line emission than 
to the radio axis. The best fit to the slope of the horizontal 
centroid has a PA of ~54°, part way between the NLR and 
ENLR axes. For the data to fit the radio axis equally well, 
either end of the outside (~25%) contour would have to be 
displaced by ~0.8". 

Fig. [j] compares the [Fe II] map direc tly with a narrow - 
band HST image of [O III] line emission ( Kaiser et al. 200C ) 
and a c ontour map of the radio jet at 8 MHz (Mundell et 



al. 1995). It can be seen that the region over which [Fell] is 



detected in the SMIRFS data corresponds to the extent of 
the narrow line region. If there were a strong correspondence 
with the radio jet, one might expect a more elongated dis- 
tribution, resolving some knot structure at 77°, where the 
outflow interacts with gas. Although detection of the larger- 
scale ENLR would have provided a better separation, the 
spatial association of [Fell] with photoionized gas is more 
convincing than with shocking by radio plasma. The F502N 
filter used for i maging [O III] covers almost all of the N LR 
velocity range (Hutchings et al. 1999; Kaiser et al. 2000), so 
represents the total emission well. It is perhaps interesting 
to note that [Fell] appears better correlated with t he po- 
sitions of low velocity dispersion clouds, plotted in ( |Kaiser 



et al. 2000), than the NLR as a whole. However, this is at 



variance with our previous suggestion that the greater ve- 
locity width of the [Fell] line compared with Pa/3 suggests 
an association with the high dispersion component. 

The Pa/3 narrow line is slightly extended in both di- 
rections, and has a faint wing out to ~1.5" South-West of 
the centre, where the NLR is brightest. Fig. ^ shows both 
the Pa/3 flux and [FeII]/Pa/3 ratio. Since Pa/3 falls off fairly 
rapidly from the centre, the ratio is highly sensitive to noise 
artifacts. The datacube was therefore smoothed with a 0.9" 
wide Gaussian kernel before producing the excitation map, 
preserving minimal resolution of spatial structure. The line 
ratio is ~0.7 at the centre and > 2 along the radio jet at 
the edges. A dip near the centre is to be expected for pho- 
toionization, since the strong incident flu x creates large full y 
ionized zones, favouring Pa/3 over [Fell] (Mouri et al. 2000). 
The high areas could reflect an enhancement of [Fell] emis- 
sion by shocks, but the line images are more indicative of a 
drop in Pa/3 at those points, with respect to [O III] , than of 
strong [Fe II] . To the South-West, a few highly-ionized clouds 
may just happen to give stronger Pa/3 away from the radio 
axis, whilst to the North-East, the NLR is angled slightly 
more towards the radio jet anyway. 



6 CONCLUSIONS 

The distribution of [Fell] emission at the centre of 
NGC 4151 has been mapped in two dimensions, using the 
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Figure 6. Map of [Fell] flux, with measurements of sample spectra. The axes marked are from Winge et al. (1999) for the ENLR and 
Kaiser et al. (2000) for the NLR and radio emission. 
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Figure 7. (left) HST image of NGC 4151 in [O III] emission, overlaid with radio contours, from Kaiser et al. (2000). (right) SMIRFS-IFU 
image of [Fe II] on the same scale. Dotted lines mark the [Fe II] axis. 
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Figure 8. (left) Pa/3 emission 

innovative SMIRFS-IFU. Although the line is only detected 
over a few spatial resolution elements, it is evident that its 
image profile is better aligned with the visible narrow line 
region than with the radio jet. The study therefore suggests 
that [Fell] arises primarily through photoionization of gas, 
by collimated X-rays from the Seyfert nucleus. The mean ve- 
locity field is consistent with this interpretation, indicating 
that [Fell], Pa/3 and the NLR all have similar bulk kine- 
matics. Likewise, the variation in [Fell] /Pa/3 is compatible 
with photoionization dominating. As suggested by previous 
authors, shock excitation by outflowing radio plasma may 
make a secondary contribution to the emission. However, 
the evidence for this comes mainly from line widths, whereas 
tentative measurements suggest, if anything, that [Fell] is 
broader along the NLR than the radio jet. 

Despite the small field and prototypical nature of the 
IFU, these results demonstrate the utility of integral field 
spectroscopy in investigating the properties of active galax- 
ies. Without a two dimensional field, it would have been 
very difficult to determine the spatial orientation of the line 
emission. Narrow-band imaging can also fulfil this require- 
ment, given appropriate filters, but follow-up spectroscopy 
is then required to form a complete picture. The experience 
gained with the SMfRFS-fFU and associated data analysis 
has paved the way for subsequent projects, both technical 
and observational. 

Reduction of the fFU data has highlighted a few impor- 
tant observational issues. In particular, location of spectra 
on the detector is difficult, and it is recommended that more 
than one reference gap between fibres should be available in 
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(right) Ratio of [Fell] /Pa/3. 

future work. With a number of sources of error and inho- 
mogeneity, the availability of multiple observations at each 
position has been very valuable. 

Whilst the results presented here support the argument 
for [Fe II] originating in photoionized gas, there is no doubt 
that observations at high spatial resolution, as well as of the 
wider ENLR, will provide further insight into the nature of 
NGC 4151 and other Seyfert galaxies. With a bright, point- 
like nucleus, NGC 4151 is an ideal target for J-band obser- 
vations with adaptive optics. This would provide a useful 
counterpart to the Hubble Space Telescope observations of 
[O III] emission, which have revealed the cloud structure and 
detailed kinematics of the NLR. In conjunction with recent 
models of photoionized and shock excited gas, comparison 
of optical and infrared lines in individual clouds would pro- 
vide important information about the physical conditions 
and kinematic structure of the NLR. The misalignment with 
the radio jet will again help determine its influence on the 
ambient gas. 
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